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THE  THERMAL  RESPONSE  OF  GYPSUM-PANEL/STEEL-STUD  WALL  SYSTEMS 
EXPOSED  TO  FIRE  ENVIRONMENTS  - A SIMULATION  FOR  USE  IN  ZONE-TYPE 

FIRE  MODELS 

Leonard  Y.  Cooper 


ABSTRACT 

This  work  develops  a method  for  simulating  the  thermal  response  of  fire-environment-exposed  wall 
systems  constructed  of  arbitrary-thickness  gypsum  panels  mounted  on  either  side  of  vertical  steel 
studs.  The  studs,  separated  at  regular  intervals,  form  an  unfilled  air  gap  between  the  panels 

The  main  objective  is  an  experimentally-validated,  modular,  thermal-wall-model  algorithm  that  can 
be  easily  integrated  into  zone-type  compartment  fire  models  and  that  can  be  used  in  “stand-alone” 
analyses. 

The  algorithm  solves  an  initial-value/boundary-value  problem  for  the  temperature  distribution 
through  the  thicknessess  of  the  two  panels  and  within  a specified  time  interval.  The  analysis  is  based 
on  temperature-dependent  thermal  properties  for  the  gypsum.  The  initial  temperature  distribution 
and  the  type  of  boundary  conditions  at  the  outer  surfaces  are  user-specified.  A variety  of  choices  for 
the  boundary  conditions  are  available.  These  include  boundary  conditions  that  are  expected  to 
satisfy  the  requirements  of  any  zone-type  fire  model  and  those  that  can  be  used  to  determine  fire 
performance  of  the  wall  systems  in  ASTM  El  19  or  ISO  834  tests.  The  algorithm  output  includes 
the  final  temperature  distribution  and,  when  outer  surface  temperatures  are  specified,  the  final  rate 
of  heat  transfer  to  these  surfaces. 


Keywords:  Algorithms,  ASTM  El  19,  compartment  fires,  fire  barriers,  fire  models,  gypsum 

board,  steel  studs,  walls,  zone  models 


INTRODUCTION  - THE  WALL  SYSTEM  AND  BASIC  ELEMENTS  OF  A THERMAL 
MODEL 


The  purpose  of  this  work  is  to  develop  a method  for  simulating  the  thermal  response  of  wall  systems 
constructed  of  gypsum  panels  mounted  on  steel  studs  to  fire  environment  exposures.  The  main 
objective  is  an  experimentally-verified,  modular,  thermal-wall-model  algorithm  that  can  be  easily 
integrated  into  zone-type  compartment  fire  models  and  that  can  be  used  in  “stand-alone”  analyses. 

Figure  1,  adopted  from  References  [1]  and  [2],  is  a sketch  of  example  wall  system  designs  under 
consideration.  In  general,  two  arbitrary-thickness  gypsum  wall  panels  are  mounted  one  on  either  side 
of  an  array  of  vertical  steel  studs.  In  practice,  each  of  the  two  panels  shown  can  involve  a single 
thickness  of  gypsum  board  or  a sandwich-type  multiple-thickness  design  of  two  or  more  well- 
contacted  boards.  Figure  1 illustrates  two  possible  assembly  designs.  One  of  these  is  referred  to 
as  a lxl-t\pe  assembly,  since  each  of  the  two  panels  involves  a single  layer  of  gypsum  board.  The 
other  is  a 1x2  assembly,  since  one  panel  is  a single  layer  of  gypsum  board  and  the  other  involves  a 
two-layer  construction.  The  studs,  separated  at  regular  intervals,  form  an  unfilled  air  gap  between 
the  panels.  As  is  the  case  in  practical  implementations  of  these  kinds  of  wall  systems,  the  spacing 
of  the  studs  is  several  times  the  thickness  of  the  air  gap.  Also,  the  studs  are  typically  fabricated  from 
relatively  thin-gage  steel  (the  studs  used  in  the  experimental  study  of  [1]  were  0.46  mm  thick)  and 
they  are  not  effective  as  paths  for  conductive  heat  transfer  between  the  panels. 

Extensive  thermocouple  data  on  the  thermal  response  of  a Figure- 1 lx2-type  wall  system  to  ASTM 
El  19  [3]  standard-fire  furnace  exposures  were  acquired  and  presented  in  Reference  [1].  These  data 
indicate  that  the  temperature  distribution  in  the  gypsum  panels,  even  relatively  close  to  the  steel 
studs,  were  substantially  one-dimensional  through  the  thickness  of  the  panels.  Consistent  with  this 
finding  and  with  an  experimentally-validated  thermal  response  model  presented  in  [2]  (all  gypsum 
panels  used  in  [1]  and  [2]  were  of  fire-rated  material),  it  will  be  generally  assumed  here  that:  1) 
relative  to  effects  of  conduction  heat  transfer,  the  steel  studs  simply  act  as  thermally  insulating 
spacers  for  the  gypsum  panels,  and  2)  radiation  exchanges  across  the  air  gap,  between  the  facing 
surfaces  of  the  gypsum  panels,  can  be  well-predicted  by  an  analysis  involving  radiative  exchange 
between  two  infinite  parallel  planes,  i.e.,  the  steel  studs  do  not  have  a significant  effect  on  modifying 
the  radiation  exchange  between  the  facing  panel  surfaces.  In  the  analysis  to  be  presented  here,  the 
time-dependent  thermal  response  of  the  gypsum  panels  will  be  simulated  by  an  idealized  system 
involving  two  initially  uniform-temperature  vertical  gypsum  board  panels,  infinite  in  extent  and 
separated  by  an  air  gap,  where  the  system  is  always  heated  at  the  two  bounding  outer  surfaces  by 
spatially-uniform  heat  fluxes. 


THE  GEOMETRY  OF  THE  MODEL  WALL  SYSTEM 

Figure  2 defines  the  generic  geometry  of  the  two-panel  wall  system  under  consideration.  As  seen 
in  the  figure,  x is  the  distance  measured  normal  to  the  panels,  from  left  to  right.  The  left-  and  right- 
hand  panels  are  designated  as  panels  1 and  2,  respectively,  and  x = x,  = 0 is  taken  as  the  left  surface 
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of  panel  1,  the  other  three  surfaces  being  located  at  x = X2,  X3,  and  X4,  where  Xj  = 0 < X2  < X3  < 
X4.  The  thicknesses  of  panel  1,  the  air  gap,  and  panel  2,  designated  as  Lj,  L^ap,  and  L2,  respectively, 
are: 


L,  = X,  - X, 


Lgap  ~ 


X,  - Xn 


L2  = 


X,  - X, 


(1) 


The  enclosed  rooms  or  open  spaces  to  the  left  and  right  of  panels  1 and  2 are  designated  as  spaces 
1 and  2 respectively.  The  surface  located  at  x = x^  is  designated  as  surface  i. 


HEAT  TRANSFER  TO  THE  PANEL  SURFACES 
Define  q"  as  the  flux  of  heat  transfer  to  a panel  at  surface  i. 

Outside-Facing  Surfaces  1 and  4;  General  Considerations 

The  two-panel  wall  system  being  considered  can  be  used  as  a barrier  that  separates  enclosed  spaces 
of  a multi-room  facility,  or  as  a barrier  that  separates  an  inside  space  and  the  outside  environment. 
Either,  or  both  outward-facing  surfaces  of  the  wall  system,  surfaces  1 or  4,  can  be  exposed  to  a fire- 
generated environment  or  to  an  ambient  environment.  In  general,  q",  i = 1 or  4,  will  be  a function 
of  Tj  = T(x  = Xj,  t),  Cj,  the  effective  gray-body  emissivity  of  surface  i,  and  of  the  variables 
that  describe  the  environment  in  the  space  to  which  surface  i is  exposed. 


q • = FjCTj,  8^;  variables  of  the  environment  to  which  surface  i is  exposed),  i = 1 or  4 (2) 


When  spaces  1 and/or  2 involve  fire  environments,  the  specification  for  Fj  and  F4  will  depend  on 
the  method  chosen  to  describe  or  simulate  the  fire  environment  there,  e.g.,  a method  consistent  with 
a particular  room  fire  model  adopted.  To  use  the  thermal  wall  model  in  a fire  model  simulation,  it 
will  be  necessary  to  provide  explicit  functional  relationships  for  the  Fj. 

Outside-Facing  Surfaces  1 and  4;  Simulating  Furnace  Test  Exposures 

This  work  describes  a thermal  model  for  Figure  1-type  wall  systems  and  verifies  it  using 
experimental  data  from  ASTM  El  19  standard  fire  exposures  as  presented  in  [1]  and  [2].  For  the 
verification,  it  will  be  necessary  to  provide  relationships  for  F^,  i = 1 and  4 of  Eq.  (2),  when  the 
environment  in  one  of  these  spaces,  say  space  1,  corresponds  to  that  found  in  a wall-furnace 
enclosure  during  a standard  fire  exposure,  and  when  the  environment  in  the  other  space,  space  2, 
corresponds  to  that  found  in  the  laboratory  space  containing  the  wall  furnace. 
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Heat  transfer  to  the  furnace-fire-exposed  surface.  Consider  the  rate  of  heat  transfer  to  the 
fumace-fire -exposed  surfaces  of  construction  elements  undergoing  ASTM  El  19-type  standard  tests, 
where  a\  erage  furnace  gas  temperatures  are  specified  to  follow  the  standard  time-temperature  curve, 
Tsto(0,  ^ specified,  for  example,  by  ASTM  El  19  [3]  or  ISO  834  [4].  It  is  well-accepted  that  Fj, 
corresponding  to  such  furnace  test  exposures,  can  be  simulated  by  (see,  e.g.,  [5]  or  [6]) 


Qi  — El  — — hp|jj^(TsjQ  - T,)  + [ci/(l/8pui^  + 1/8]  - 1)](T^td  - T{) 


(3) 


where  8pu]^  is  a combined,  effective,  gray-body,  surface  emissivity  of  the  furnace  walls  and  fire  gases, 
8]  is  the  effective  gray-body  emissivity  of  surface  1,  and  o,  the  Stephan-Bolzmann  constant,  and 
hp^TR,  the  heat  transfer  coefficient  for  convective  transfer  to  the  furnace-exposed  surface  of  the  system 
being  tested,  are  given  by 


= 25  W/(m-K);  o = 5.7(10-*)  W/(m^K") 


(4) 


On  the  right  side  of  Eq.  (3),  the  first  term  provides  a rough  estimate  of  the  convective  component 
of  the  heat  transfer.  It  is  a result  of  the  relatively  vigorous  forced  convection  flows  which  are 
characteristic  of  a furnace  environment.  The  second  term  is  the  radiative  component  of  the  heat 
transfer.  A precise  estimate  of  the  convective  heat  transfer  is  not  required  since,  for  furnace 
temperatures  and  (TpL-p,  - T,)  temperature  differences  of  interest,  the  radiation  term  is  typically  much 
larger  than  convection  term. 

Heat  transfer  to  the  laboratory-exposed  surface.  Heat  transfer  to  the  ambient-environment- 
exposed  surfaces  would  also  have  convective  and  radiative  components.  The  latter  would  be 
characterized  by  radiative  exchanges  between  surface  4 and  a black-body  far  field  radiating  at  the 
ambient  temperature  T^g-  Assuming  a relatively  quiescent  laboratory  environment,  any  convective 
heat  transfer  to  surface  4 would  be  a result  of  natural  convection. 

.A.n  estimate  for  heat  transfer  to  a vertical  surface  of  height  H at  uniform  temperature  T4  from  a 
quiescent  atmosphere  at  T^^ib  < ^4  t>e  obtained  with  the  use  of  the  heat  transfer  coefficient  [7] 


h^„3=  1.3(T,-T^,B)‘'’W/(m'K^'') 


(5) 


for  a turbulent  boundary  layer  flow,  provided 


Ti  - > 50  K/(H/m)^ 


(6) 
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Eq.  (5)  will  be  used  below.  At  times  when  Eq.  (6)  is  not  satisfied,  it  will  be  assumed  that  both  acmal 
and  (poorly)  estimated  convective  heat  transfer  at  surface  4 would  have  a negligible  effect  on  the 
overall  thermal  response  of  the  wall  system. 

In  view  of  the  above,  under  standard  fire  test  conditions  F4  will  be  taken  to  be 


Q4  - F4  - " hAMBCf'4  ■ Tamb)  ■ ^4  ^(^4  - T^b) 


(7) 


Inside-Facing  Surfaces  2 and  3;  Simulating  Heat  Transfer  Across  the  Air  Gap 

As  differences  between  T2  and  T3  begin  to  develop  during  the  course  of  a fire  exposure,  q^'  and  q^’ 
will  have  a radiative  component,  q'^  ^nd  q's  RAO’  occurring  as  a result  of  radiation  ex- 
changes through  the  air  layer,  and  a component,  q2,coN  Qs.con’  resulting  from  surface-to-air 
conduction/convection  exchanges. 


Q2  ~ Q2,RAD  Q2,C0N’  Q3  ~ ^3.RAD  + ^3, CON 


(8) 


For  the  purpose  of  calculating  the  radiation  exchanges  it  is  assumed  that  the  air  layer  is  transparent. 
Therefore  [8] 


42,RAD  ~ [o/(1/82  + 1/^3  - 1)](T^  - — - Q3,raD 


(9) 


Assume  that  the  heat  capacity  of  the  mass  of  air  in  the  gap  is  negligible  in  the  sense  that  a relatively 
small  rate  of  energy  transfer  is  required  to  maintain  the  average  air  temperature,  Tair,  at 


Tair  = (Ti  + T2)/2 


(10) 


Then,  in  terms  of  an  effective  thermal  conductivity,  k^pp,  the  conduction/convection  components  can 
be  estimated  by  the  simple  cross-gap-type  conduction  model  [8] 


Q2,C0N  ~ (1^EFf/Lgap)(T3  " T2)  - - q3  c 


(11) 
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where 


^EFF^.AjR  “ f(GrPr);  Gr  — g(T3  - T2)L(3^(T^v  );  Pr  ~ 0.7  ■ (12) 


In  Eq.  (12),  Gr  is  a Grashoff  number  based  on  the  air  properties  and  the  air  gap  thickness,  g is  the 
acceleration  of  gravity,  and  Pr  and  v are  the  Prandtl  number  and  kinematic  viscosity,  respectively, 
of  air  at  temperamre  T^.  f(GrPr),  from  a correlation  of  experimental  data  [9],  is  plotted  in  Figure 
1 1-14  of  [8].  A curve  fit  for  these  data  is  [10] 


1, 


log(GrPr)  <3; 


f(GrPr)  = 


O.nCGrPrf”, 


3.8  < log(GrPr)  < 6; 


0.40(GrPr)°  “°,  6 < log(GrPr)  < 8 


(13) 


f(GrPr)  of  Eq.  (12)  is  such  that  f is  monotonically  increasing  with  GrPr,  approaching  15.8  at  GrPr 
= 10^.  For  problems  of  interest  here,  GrPr  will  likely  never  exceed  10^,  corresponding  to  f = 10,  and 
k^  would  fall  in  the  range,  0.035  W/(mK)  < k^  < 0.070W/(mK) . For  example,  using  Eqs.  (10)- 
(13)  consider  test  data  presented  in  Figure  11  of  [1]  (L^^  = 0.09  m)  at  the  relatively  early  and  late 
times,  10  and  30  min,  respectively: 


at  t = 10  min:  T.  = 380  K;  T3  = 355  K: 

- 368  K;  k^R=  0.031  W/(mK);  3.2(10-^)  mVs; 

GrPr  = 3.3(10^);  f = = 4.4;  q”  ^on  = 38  W/m^  q2.RAD  = 280  W/m' 

(14) 

at  t = 30  min:  = 790  K;  T3  = 690  K: 

T^  = 740K;  k^  = 0.054  W/(mK);  7.2(10-^)  mVs; 

GrPr  = 1.3(10^),  f = k^pp/k^ij^  = 3.3;  q2,coN  = 198  WW;  q2.RAD  = 9300  W/m^ 


The  above  estimates  suggest  that  for  the  type  of  problem  considered  here  (i.e.,  for  wall  heating 
problems  where  temperature  increases  of  the  gypsum  panels  are  considered  to  be  significant  when 
they  are  of  the  order  of  several  hundred  Kj,  when  the  net  rate  of  heat  transfer  across  the  gap  starts 
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to  become  significant,  q'^  RAD  » q2,coN-  Based  on  this  it  is  convenient  to  use  the  approximation: 


^2  ~ “ Qs  “ 42,RAD 

A MODEL  FOR  THE  THERMAL  RESPONSE  OF  THE  GYPSUM  PANELS 
Heating  of  the  Gypsum 

At  ambient  temperature  and  prior  to  any  heating,  gypsum  panels  are  made  up  of  calcium  sulfate 
dihydrate,  CaS04-2H20,  that  contains  about  21%  by  weight  of  chemically-combined  water,  and 
additional,  relatively-small  amounts  of  absorbed  free  water  [11].  Also,  practical  Figure  1-type  wall 
systems  include  a relatively-thin  paper  covering  on  the  exposed  surfaces  of  the  gypsum  panels.  In 
the  fire  tests  of  [1],  this  paper,  approximately  0.4  mm  thick,  was  observed  to  bum  off  within  3 min 
of  the  start  of  the  test.  The  analysis  for  the  thermal  response  of  the  gypsum  panels,  to  be  developed 
here,  will  ignore  the  effects  of  such  surface-mounted  paper. 

When  the  temperature  of  a panel  is  raised  above  80  °C,  the  chemically-bound  water  begins  to  be 
released  in  a process  that  involves  two  stages.  Most  of  the  water  is  released  in  the  first  stage,  as  the 
CaS04-2H20  becomes  calcium  sulfate  hemihydrate,  CaS04-l/2H20  (plaster  of  Paris).  This  stage  is 
usually  complete  by  the  time  the  material  reaches  125  °C.  As  heating  is  continued,  the  rest  of  the 
water  is  released  and  the  remaining  material  goes  to  anyhydrous  calcium  sulphate,  CaS04.  During 
these  processes,  when  the  material  is  at  and  above  the  boiling  temperature,  absorbed  and  released 
water  goes  to  the  gaseous  phase  in  a process  that  requires  the  absorption  of  relatively  large  amounts 
of  energy.  The  reader  is  referred  to  Reference  [1 1]  for  more  detail  of  these  processes,  and  for  further 
reference  to  other  relevant  literature. 

When  the  steam  is  released  in  depth,  moisture  migration  must  occur  before  it  can  actually  reach  one 
or  the  other  bounding  surfaces  and  be  released  to  the  adjacent  atmosphere.  Also,  the  migration  of 
moisture  toward  the  relatively-cool  in-depth  portion  of  the  gypsum  can  lead  to  recondensation  of  the 
water  vapor,  which  will,  in  turn,  lead  to  some  modification  locally  to  the  permeability  of  the 
gypsum.  Some  details  of  modeling  the  moisture  migration  are  presented  in  Reference  [12]. 

Modeling  the  Thermal  Response  of  the  Gypsum  and  Its  Thermal  Properties 

Both  gypsum  panels  of  a particular  wall  system  design  are  assumed  to  have  identical  properties. 

A conduction  heat  transfer  analysis  of  the  thermal  response  of  Figure  1-type  wall  systems  to  ASTM 
El  19-type  furnace  heating  is  presented  in  [2].  The  analysis  uses  a temperature-dependent,  specific- 
heat  model  for  the  gyspum  panels  that  is  deduced  directly  from  differential  scanning  calorimeter 
measurements,  at  2°C/min,  on  a relatively-small  40  mg  sample.  Note  that  representation  of  specific 
heat  in  the  conduction  equation  based  on  such  small-sample  data  is  equivalent  to  an  assumption  that 
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1)  the  thickness  of  the  panels  are  small  enough  for  substantially  unrestricted  moisture  migration  to 
occur  and  2)  any  effects  of  recondensing  water  vapor  are  not  significant. 

As  reponed  in  [2],  and  especially  at  times  w'hen  T2  and  T3  exceed  approximately  200  °C,  wall  model 
simulations  of  ASTM  El  19-type  fire  resistance  tests  were  shown  to  compare  favorably  with 
expenmental  data.  Expenmental  data  presented  in  [2]  include  the  above-mentioned  test  data  of  [1]. 
.All  these  data  were  acquired  with  g>psum  panels  and  gypsum  samples  of  fire-rated  Type  X material. 

Reference  [11]  develops  a two-dimensional  heat  conduction  model  for  a Eigure  1-type  wall  system 
(same  thickness  panels  as  in  [2]),  w'here  wood  studs  replace  steel  studs  and  where  the  thermal 
response  of  the  wood  studs  is  an  important  and  desired  feature  of  that  analysis.  In  the  part  of  the 
model  that  deals  with  the  response  of  the  gypsum  panels,  it  is  assumed,  in  contrast  to  [2],  that 
significant  moisture  migration  does  not  occur.  In  particular,  subsequent  to  the  phase  change  of 
released  water  from  liquid  to  steam,  Cp(T)  is  determined  from  an  analysis  that  treats  the  gypsum  as 
a combined  solid/steam  system.  Einally,  in  [11],  p(T)  is  determined  from  thermogravimetric 
analyzer  data  acquired  at  a scanning  rate  of  20  °C/min  and  k(T)  from  steady-state  measurements. 

As  reported  in  [13],  and,  again,  especially  at  times  when  T2  and  T3  exceed  approximately  200 
°C,  simulation  results  using  an  advanced  version  of  the  Reference-[1 1]  model  were  shown  to 
compare  favorably  with  experimental  data  from  ASTM  El  19-type  fire  resistance  tests. 

Wall  models  like  those  of  [2],  [1 1],  and  [13]  that  depend  on  a straight-forward  conduction  equation 
analysis  and  use  of  effective,  temperature-dependent  thermal  properties  can  only  be  expected  to  yield 
accurate  results  to  the  extent  that  details  of  moisture  migration  phenomena,  including  the  details  of 
local  moisture-dependent  property  variations  are  not  significant.  The  favorable  results  mentioned 
above  of  the  models  of  [2],  [1 1],  and  [13],  suggest  that  when  the  gypsum  exceeds  approximately  200 
°C,  moisture-dependence  of  the  properties  and  moisture  migration  details  can  be  so  ignored. 

The  reader  is  referred  to  References  [14],  [15],  and  [16]  for  additional  insight  on  modeling  the 
thermal  response  of  the  gypsum  and  its  thermal  properties. 

Representations  for  the  Density,  Specific  Heat,  and  Thermal  Conductivity  of  Gypsum 

In  modeling  the  thermal  response  of  Eigure- 1 -type  wall  systems  to  fire  exposures,  especially  in 
applications  of  present  interest  where  simulations  to  temperatures  of  200  °C  and  higher  are  of 
particular  interest,  it  will  be  assumed  here,  consisent  with  the  above  discussion,  that:  1)  the  thermal 
response  of  the  g>'psum  panels  can  be  deduced  from  a traditional  heat  conduction  model  that  uses 
p(T),  Cp(T),  and  k(T)  representations  determined  with  the  use  of  traditional  property-measurement 
techniques,  and  2)  the  details  of  moisture  migration  can  be  ignored. 

The  wall  model  to  be  developed  here  will  be  verified  with  the  wall  furnace  test  data  of  [1]  and  [2]. 
p(T),  Cp(T),  and  k(T)  of  the  gypsum  panels  used  in  those  tests  were  measured  under  quasi-steady- 
state  temperature  conditions  and  reported  in  [2]  as  linear  interpolations  between  [property. 
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temperature]  data  pairs. 

Property  representations  used  here  for  wall-model  validation  and  presented  in  APPENDIX  A are 
based  on  the  property  representations  of  [2].  The  p(T)/p(T  = 20  °C),  Cp(T),  and  k(T)  representations 
are  plotted  in  Figures  3-5,  respectively. 


THE  INITIAL/BOUNDARY-VALUE  PROBLEM  FOR  THE  THERMAL  RESPONSE  OF 
THE  WALL  SYSTEM 

The  objective  of  this  work  is  an  algorithm  for  determining  the  temperature  distribution  T(x,  t)  in  the 
gypsum  wall  panels  of  a Figure-2-type  wall  system  for  any  time  t,  tj  < t < t2, 
for  specified  initial  condition,  Tp^Cx)  = T(x,  tj).  The  algorithm  must  solve  the  following  generic 
initial/boundary-value  problem: 


Find  T(x,  t)  for  tj  < t < t2  and  Xj  <x  < X2  or  X2  < X3  < x < X4  satisfying 

p(T)Cp(T)aT/at  = k(T)a'T/ax2  + (ak/ax)(aT/ax)  ( 1 6) 

subject  to  the  initial  conditions: 


T(x,  t — tj)  — Tj]^(x)  (17) 

for  specified  Tn,^(x),  and  the  boundary  condition: 

at  X = Xi:  T = T,(t)  or  k(T)aT/ax  = - q';(t)  (18) 


k(T)aT/ax 


= k(T)aT/ax 


= [o/(1/82  -h  1/03-  1)][T(X3,  t)^  - T(X2,  t)1  (19) 


X = X2 


X = X3 


at  X = X4:  T = T4(t)  or  k(T)aT/ax  = q:;(t) 
for  specified  Ti(t)  or  q'i'(t)  and  T4(t)  or  q4(t). 


(20) 


In  the  above,  note  that  Eq.  (19)  follows  from  Eqs.  (9)  and  (15). 

When  the  specific  problem  of  interest  is  to  find  the  response  of  the  wall  system  to  a standard  fire 
furnace  exposure,  the  problem  involves  the  following  explicit  conditions: 


When  surface  1 is  exposed  to  a standard  fire  furnace  exposure,  the  initial  condition  of  Eq. 
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(17)  becomes 


T(x,  t — 0)  — T^ib  — T5yQ(t  — 0)  (17) 

for  specified  T^mb’  ^1^^  boundary  conditions  of  Eqs.  (18)  and  (20)  become 

at  X = X,:  k(T)cT/ax  = - hp^,,R(TsTn  - T)  - [o/(l/£p„R  + 1/e,  - DJCTstd  ' T^)  d^') 

at  X = X4;  k(T)cT/0x  = - h^.^^B(T  - T^g)  - 840(1^  - T^g)  (20 ) 

where  recommended  values  for  and  h^^g  are  given  in  Eqs.  (4)  and  (5),  respectively. 

THE  INTTIAI7BOUNDARY  VALUE  PROBLEM  IN  TRANSFORMED  VARIABLES 

The  algorithm  for  solving  Eqs.  (16)-(20)  will  use  a computer  software  solver  that  motivates 
redefinition  of  the  problem  in  terms  of  the  following  transformed  spatial  variables  and  new 
dependent  and  independent  variables: 


For  X,  <x  < X2: 


j(i)  ^ j.  0x(i)/ax  = ai/ax;  X = X - X,  (i.e.,  0 < X < L,); 
X,  = = l/[p(r^))Cp(T(^^)];  Kj=  K,(r‘))  = 

For  X3  < X < X4: 


(21) 


= T,  ar^vax  = - (L,/L,)aT/ax;  X = LJI  - (x  - x^)/L,]  (i.e.,  0 < X < L,); 

(22) 

A..  = = (L,/L0V[p(r'^)Cp(r'^)];  K.  = 


In  the  transformed  variables  of  Eq.  (22):  the  original  depth  of  panel  2,  L2,  is  expanded/contracted  to 
L[,  the  depth  of  panel  1;  the  laboratory-exposed  surface  of  panel  2,  originally  at  its  right  side,  is 
transformed  to  its  left  side,  at  the  new  origin,  X = 0;  and  the  cavity-exposed  surface  of  panel  2, 
originally  on  the  left,  is  transformed  to  the  right  side,  at  X = Lj.  The  transformation  of  Eq.  (22)  is 
depicted  in  Figure  2. 

Then  the  problem  of  Eqs  (16)  - (20)  becomes: 


10 


Find  T‘'’(X,  t)  and  T®(X,t)  for  t,  < t < and  0 < X < L,  satisfying 
aT<‘Vat  = K,(T<‘>)a^T“Vax^  + A,|(T<‘>)(ak/ax)(aT*'Vax) 

(23) 

ax'^vat  = K2(T'"')a^T<^vax2  + X2(T'^')(ak/ax)(aT<^vax) 

subject  to  the  initial  conditions 

T<‘>(X,  t = t,)  = Tj^(x  = X + X,);  T®(X,  t = t,)  = T^(x  = x,  + Lj  - XL^/L,,  t,) 

(24) 

and  the  boundary  conditions 
at  X = 0: 

= Ti(t)  or  k(r ‘^)aT<*VaX  = - q';(t);  (25) 

= T4(t)  or  k(T^'^)aT('VaX  = - q'^t)  (26) 

atX  = L,: 

k(T(n)aT^^vax  = - (L2/Li)k(r'^)aT('vax  = [a/(i/82  + 1/83-  i)](r'^  - r^^") 

(27) 

When  surface  1 is  exposed  to  a standard  fire  furnace  exposure,  the  initial  conditions  of  Eq. 
(24)  become 

t = 0)  = r2>(X,  t = 0)  = (24') 

and  the  boundary  conditions,  Eqs  (25)  and  (26),  become 

at  X = 0: 

k(T”»)aT'‘VaX  = - hpuRdsTO  - T'‘>)  - [a/d/ep,;^  + 1/e,-  - T<"")  (25') 

k(T<^>)aT«Vax  = - T^b)  + 64  a(T<^'^  - T^mb)  (26') 


SOLVING  THE  PROBLEM 


The  Subroutine  GY  PST 

The  FORTR.A.N  77  subroutine  GYPST  and  associated  computer  software  was  developed  to  solve 
the  initial/boundar>  -value  problem  of  Eqs.  (16)  - (20).  The  subroutine  uses  the  software  MOLID 
[17],  which  solves  for  T^\X,  t)  and  E‘\X,  t)  of  the  transformed  problem  of  Eqs.  (23)  - (27).  Once 
these  have  been  obtained,  then,  following  Eqs.  (21)  and  (22),  the  T(x,  t)  solution  of  the  original 
problem  is  retrieved  according  to: 


For  0 < X < L,: 


T = aT/ax  = ar ‘Vax,  x = x,  + x , i.e.,  x,  <x  < X2 

T = V-\  aT/ax  - (L2/L,)aT^‘VaX,  x = X3  + L.  - X U/L,  i.e.,  X3  <x  < X4 


(28) 


Solving  the  Transformed  Problem  with  MOLID 

The  transformed  problem  of  Eqs.  (23)  - (27)  is  solved  with  MOLID  [17].  This  is  a set  of 
FORTRAN  subroutines  that  provides  for  the  method-of-lines  solution  for  systems  of  initial- 
boundaiy-value  partial  differential  equations  in  one  space  dimension,  including  systems  that  involve 
linear  and  non-linear  parabolic  equations  subject  to  a wide  variety  of  boundary  conditions. 

In  MOLID,  spatial  derivatives  are  approximated  using  either  2nd,  3rd,  or  6th-order  accurate  finite 
differences  on  an  equi-spaced  grid.  (In  the  present  use  of  MOLID,  2nd-order  accurate  finite 
differences  are  used.)  The  resulting  ordinary  initial  value  problem  is  then  solved  using  the  adaptive 
integration  routine  of  [18].  In  the  calculation,  estimates  of  the  time  step  error  at  the  spacial  grid 
points  are  kept  less  than  a user-specified,  time-integrator,  local-error,  tolerance  parameter  EPSCAL. 
(In  most  example  calculations  to  be  presented  below,  EPSCAL  is  specified  to  be  0.00001.) 

Using  GYPST 

To  solve  a particular  wall-system  fire-exposure  problem  with  GYPST,  one  must  first  specify 
T(x,  t = TINT j at  each  of  NPTS  equi-distant  nodes  in  panel  1 and  NPTS  equi-distant  nodes  in  panel 
2.  Then,  a call  to  the  subroutine  GYPST  with  this  and  other  input  data  leads  to  a calculation  of  the 
thermal  response  of  the  wall  system,  subjected  to  one  of  nine  possible  sets  of  boundary  conditions, 
between  t = TENT  and  specified  time  t = LEAST.  The  set  of  boundary  conditions  of  interest  is 
identified  by  the  input  value  EBCTYP  =1-9.  The  GYPST  output  yields  the  newly  calculated 
T(x,  t = TLAST)  at  each  of  the  2NPTS  node  points  of  the  two  panels.  If  the  specified  EBCTYP  set 
of  boundary'  conditions  is  one  that  does  not  explicitly  specify  heat  transfer  rates  at  surfaces  1 and  4, 
then  the  GYPST  output  also  yields  newly  calculated  values  of  these  heat  transfer  rates. 
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The  choices  of  boundary  conditions  include  IBCTYP  = 5,  which  is  used  in  all  GYPST  simulations 
to  be  discussed  below.  This  corresponds  a simulation  of  the  scenario  where  a specified  wall  system 
is  subjected  to  an  ASTM  El  19  fire  furnace  test.  In  particular,  all  GYPST  simulations  to  be  presented 
below  will  involve  use  of  Eqs.  (17')  and  (18')  with  TgjD  calculated  according  to  [2]  , 


Tstd(0  - TastME119(0 

= Tamb  + 750.{  1.-  exp[-  3.79553(t/3600)^'2]}  + 170.41(t/3600)''^  t in  [s],  T’s  in  °C 

A concise  description  of  GYPST,  extracted  from  the  GYPST  source  code,  is  presented  in 
APPENDIX  B.  This  includes  brief  explanations  of  all  options  for  boundary  conditions  and 
definitions  of  all  parameters  of  the  subroutine  call  sequence. 


SIMULATING  FURNACE  EXPOSURE  TESTS  OF  REFERENCES  [1]  AND  [2]; 
VERIFICATION  OF  THE  WALL  MODEL 

Standard  Fire  Furnace  Tests  on  Two  Different  Figure-l-Type  Wall  Systems 

Extensive  tabulated  thermocouple  data  on  the  thermal  response  of  a full-scale,  3.0  m (width)  x 3.7 
m (height).  Figure- 1 -type  wall  system  to  standard-fire  furnace  exposures  were  acquired  and 
presented  in  Reference  [1]  (refer  to  the  test  of  Assembly  F-07).  The  wall  system  tested  was  a 1x2- 
type  assembly,  where  the  single-layer  panel  was  the  one  exposed  to  the  furnace  environment.  In  the 
particular  case  of  the  F-07  wall  design,  all  panels  were  0.0127  m thick,  i.e.,  a 0.0127  m thick  gypsum 
panel  on  the  fire-exposed  side  of  the  studs  and  a 0.0254  m-thick  two-panel-sandwich  arrangement 
on  the  other.  Assembly  F-07  test  data  are  also  presented  in  Figure  9 of  [2].  This  assembly  and  test 
will  be  referred  to  below  as  the  “1x2  assembly/test.” 

In  the  1x2  test,  the  furnace  was  operated  to  reproduce  the  CANAJLC-S101-M89  [17],  standard-fire 
temperature-time  curve,  which  is  similar  to  ASTM  El  19  [2]. 

In  Figure  8 of  Reference  [2]  there  are  also  presented  plotted  test  data  of  the  thermal  response  of 
another  3.0  m x 3.7  m Figure- 1 -type  wall  system,  which  was  also  subjected  to  an  ASTM  El  19-type 
exposure.  This  assembly  was  a Ixl-type  design,  where  the  panels  used  were  each  0.0159  m thick. 
This  will  be  referred  to  below  as  the  “1x1  assembly/test.” 

Location  of  Temperature  Measurements 

During  the  tests,  temperatures  of  the  assemblies  were  measured  by  thermocouples  that  were  attached 
to  gypsum  panel  and  steel  stud  surfaces  at  different  elevations  and  at  different  lateral  positions  of  the 
assembly.  Of  interest  here  are  the  temperatures  on  the  gypsum  panel  surfaces.  Generic  locations 
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of  relevant  thermocouples,  indicated  in  Figure  1,  were: 


• hire-exposed  panel: 

surface  exposed  to  the  cavity;  (as  in  \2]  refer  to  these  as)  BI7Cav.(exp.| 
Utiexposed  panel(s): 

surface  exposed  to  the  cavity;  HI7Cav.|iinexp.| 

contact  surface  between  the  two  panels  (only  lor  1x2  assembly);  HL/FL|unexp.| 

surface  exposed  to  the  ambient/laboratory,  where  the  thermocouples  were 
unprotected;  lJnExp.|hare) 

surface  exposed  to  the  ambient/laboratory,  where  the  thermocouples  were  protected 
by  insulating  pads;  UnKxp.( under  pads).  (Since  these  thermocouples  are  protected 
from  radiative  heat  transfer  exchanges,  they  presumedly  yield  measurements  of  the 
temperatures  of  the  surfaces  to  which  they  are  attached  that  are  more  accurate  than 
those  of  the  IJnKxp.lhare]  thermocouples.) 

The  reader  is  referred  to  References  1 1 ] and  [21  for  details  of  thermocouple  locations. 

Unless  stated  otherwise,  in  the  case  of  the  1x2  te.st,  experimental  data  plotted  below  for  temperatures 
at  the  above  generic  locations  will  always  refer  to  average  temperatures  as  reported  in  Table  3 ol  1 1 1 
and  as  plotted  in  Figure  9 of  |2|.  Note  that  the  data  for  BL/FL|unexp.|,  which  is  from  Table  3 of 
[ 1 j,  does  not  appear  in  Fhgure  9 of  |2j.  When  data  from  individual  thermocouples  at  the 
HL/('av.|exp.|  locations  (thermocouples  28,  29,  32,  33,  36,  and  37)  and  the  BL/Cav.[unexp.l 
locations  (thermocouples  30,  3 1 , 34,  35,  38,  and  39  ) are  plotted  (e.g.,  in  Figure  9),  these  are  taken 
from  'fable  2 of  [ 1 ]. 

In  the  case  of  the  1x1  test,  experimental  data  to  be  presented  below  will  always  refer  to  average 
temperatures  as  plotted  in  Figure  8 of  [2j. 

From  the  above  referenced  data,  T^;^„  for  the  1 x2  and  1 x 1 test  are  taken  to  be  26  ‘’C  and  23  “ C, 
respectively. 

Simulating  the  1x1  and  1x2  lests  With  (iYPST 

In  addition  to  the  already-specified  parameters  of  the  1x1  and  1 x2  tests,  GYFST  will  require  inputs 
for  e,,  i = 1 - 3,  and  NP'fS.  'fhe  effect  of  varying  these  parameters  will  be  studied  in  the 

various  simulations  to  be  reported  below.  All  specified  parameters  of  the  1x2  and  1x1 
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assemblies/tests  and  their  GYPST  simulations  are  summarized  in  Table  1. 


Note  in  Table  1 that  in  the  1x2  assembly/test  simulations,  odd  numbers,  1 1 and  21,  were  chosen  for 
NPTS.  This  was  done  for  convenience,  so  that  for  the  unexposed  panel,  panel  2,  the  location  of  the 
middle  node  point  would  correspond  exactly  to  the  location  of  the  contact  surface,  between  its  two 
“sandwiched”  gypsum  panels,  where  the  BL/FL[unexp.]  thermocouple  data  were  acquired. 

Other  Details  of  the  Simulations 

For  every  simulation  to  be  presented,  the  simulation  test  time  was  set  at  4200  s.  The  wall-system 
reponse  for  the  entire  duration  of  the  test  was  determined  from  140  separate  30  s calls  to  GYPST 
(i.e.,  TLAST  = TINT  + 30.),  where  the  temperature  output  at  the  model  node  points  from  one  30  s 
simulation  interval  was  used  as  the  temperature  input  at  these  same  node  points  for  the  next  30  s 
interval.  All  calculations  were  carried  out  on  a PC  with  a 90  MHz  Pentium  processor. 

Expected  Limits  for  the  GYPST  Simulations 

For  both  the  1x2  and  1x1  assemblies/tests,  the  measured  temperatures  of  the  two  cavity-exposed 
surfaces,  surfaces  2 and  3 became  nearly  equal,  in  an  abrupt  manner,  as  they  reached  approximately 
600  °C.  This  was  interpreted  in  [2]  as  indication  of  a failure  in  panel  1,  and  led  to  the  conclusion 
that,  for  the  two  tests  reported,  the  “gypsum  board  is  no  longer  in  place  when  [throughout  its 
thickness]  its  temperature  exceeds  600  °C.”  Since  the  wall  model  can  not  be  expected  to  provide 
valid  simulations  when  one  of  the  panels  is  “no  longer  in  place,”  e.g.,  when  a portion  of  the  fire- 
exposed  panel  cracks  and  falls  away,  the  validity  of  the  present  GYPST  simulations  cannot  be 
expected  to  go  beyond  the  time  when  the  temperature  of  surface  2,  BL/Cav.[exp.],  reaches  600  °C. 


RESULTS  OF  GEPST  SIMULATIONS  OF  THE  1X1  ASSEMBLY/TEST 
This  section  presents  results  of  GYPST  simulations  of  the  1x1  assembly  test. 

GYPST  Simulations  Using  Different  Values  of  Emissivity 

No  guidance  is  available  on  the  most  appropriate  values  for  characterizing  the  e^’s  for  the  gypsum 
surfaces.  Similarly  unavailable  is  the  effective  value  for  the  particular  furnace  used  to  carry  out 
the  tests.  Therefore,  GYPST  simulations  of  the  test  included  a parameteric  study  involving  several 
different  sets  of  values  for  these  emissivities.  Two  series  (two  different  values  of  NPTS)  of  four 
GYPST  simulations  each  (different  values  of  and  Ej’s)  were  carried  out.  In  three  of  the  four 
latter  simulations,  and  all  q ’s  were  assumed  to  be  identical,  0.9,  0.8,  or  0.7.  In  the  fourth 
simulation,  all  €j’s  were  assumed  to  be  0.5  and  €puR  was  assumed  to  be  0.8  (since 
an  assumed  value  of  Epu^  = 0.5  is  considered  to  be  unreasonably  low).  The  two  series  involved 
different  levels  of  calculation-grid  refinement,  NPTS  = 10  or  20.  A third  series  considered  the  effect 
of  changing  the  specified  numerical  precision  parameter,  EPSCAL. 
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Simulations  with  NPTS  = 20:  Results  of  GYPST  temperature  simulations  with  NPTS  = 20  and 
the  different  values  of  emissivity  are  plotted  with  measured  temperatures  in  Figures  6a-d.  Included 
in  the  plots  are:  the  calculated  temperatures  for  the  furnace-exposed  surface  of  panel  1,  designated 
as  Exp.;  the  temperamre  of  the  furnace,  as  determined  from  Eq.  (29);  and  the  temperatures 
calculated  at  locations  corresponding  to  the  locations  BL/Cav.[exp.],  BL/Cav.[unexp.],  and 
UnExp.[under  pads]  ( or  UnExp.[bare]). 

As  can  be  seen  in  Figures  6a-c,  when  the  and  the  e^’s  are  0.9,  0.8,  or  0.7,  the  comparisons 
between  simulated  and  measured  temperature  are  ver>'  favorable  (up  to  the  time  when  the  calculation 
is  invalidated  because  of  panel  1 failure,  when  the  measured  temperature  at  BL/Cav.[exp.]  reaches 
approximately  600  °C),  with  comparisons  improving  with  increasing  emissivity.  Of  the  three 
simulations,  the  two  corresponding  to  ep^R  = = 0.9  and  0.8  (Figures  6a  and  6b)  lead  to  the  best 

comparisons  between  measured  and  calculated  values,  with  little  variation  between  the  results  in 
these  two  cases.  In  general,  the  dependence  of  simulated  temperatures  on  emissivity  is  seen  to  be 
relatively  weak  in  this  0.9  - 0.7  range  of  emissivity.  However,  a clear  degradation  in  the  comparison 
between  measured  and  calculated  temperatures  is  seen  for  the  simulation  of  Figure  6d,  where  the  Ej’s 
have  been  reduced  to  0.5. 

Calculation  times  for  all  the  above  simulations  were  approximately  160  s. 

Reducing  the  value  of  NPTS  and  EPSCAL 

If  GYPST  is  to  be  used  in  zone-type  fire  model  simulations,  then  computational  time  becomes  a 
critical  issue.  For  this  reason  it  was  of  interest  to  explore  the  effect  of  reduction  in  NPTS  and 
EPSCAL  on  simulation  performance,  in  general,  and  on  computational  time,  in  particular. 

Simulations  with  NPTS  = 10:  The  previous  simulations  were  carried  out  again  with  all  parameters, 
except  for  NPTS,  identical  to  those  used  previously.  For  this  case  NPTS  was  reduced  to  10. 

Calculation  times  for  all  the  new  simulations  were  reduced  from  approximately  160  s to 
approximately  25  s.  The  results  of  the  particular  simulation  that  corresponds  to  that  of  Figure  6a  are 
presented  in  Figure  7.  As  can  be  seen  from  these  two  figures,  the  calculated  temperatures  for  the 
new  NPTS  = 10  calculations  do  not  compare  as  favorably  to  measured  values  as  do  the  NPTS  = 20 
calculation  results.  Also,  use  of  the  relatively-coarse  NPTS  = 10  grid  introduces  “wiggles”  into  the 
calculated  temperature  histories.  Nevertheless,  at  a significant  savings  in  computational  time,  the 
new  calculation  provides  results  which  still  compare  relatively  well  with  the  measured  temperatures. 

In  the  case  of  the  Figure  6b-d-type  simulations,  the  comparison  of  the  NPTS  = 10  to  the  NPTS  = 20 
results  were  similar  to  the  Figure  6a/Figure  7 result. 

Effect  of  increasin2  EPSCAL:  Figure  6a-type  simulations  were  carried  out  with  all  parameters, 
except  for  EPSCAL,  identical  to  those  used  in  the  original  simulation.  EPSCAL  was  increased  from 
0.(XXX)1  to  0.0001,  and  then  to  0.001,  and  .01.  The  calculation  times  for  these  three  new  simulations 
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remained  in  the  range  160-155  s. 


Results  for  the  EPS  CAL  = 0.01  case  are  plotted  in  Figure  8.  As  can  be  seen,  except  for  the  predicted 
temperatures  Exp.  of  the  furnace-exposed  surface  of  panel  1 , where  scatter  of  computed  temperatures 
within  a relativly  narrow  range  can  be  observed,  the  difference  between  these  results  and  those  of 
Figure  6a,  for  EPSCAL  = 0.00001,  can  not  be  detected  visually. 

When  EPSCAL  was  increased  to  0. 1 , after  several  of  the  30  s-interval  time  steps,  the  calculation 
resulted  in  an  “overflow,”  and  the  4200  s simulation  was  not  completed  successfully. 


RESULTS  OF  GYPST  SIMULATIONS  OF  THE  1X2  ASSEMBLY/TEST 
This  section  presents  results  of  GYPST  simulations  of  the  1x2  assembly  test. 

GYPST  Simulations  Using  Different  Values  of  Emissivity 

Simulations  were  carried  out  for  the  1x2  assembly/test  with  the  same  sets  of  e’s  as  for  the  1x1 
assembly/test,  and  with  two  different  values  of  NPTS.  Here,  NPTS  was  chosen  to  be  11  or  21. 

Simulations  with  NPTS  = 21:  Results  of  GYPST  simulations  with  NPTS  = 21  and  the  different 
values  of  emissivity  (corresponding  to  those  of  Figures  6a-d)  are  plotted  with  measured  temperatures 
in  Figures  9a-d.  Included  in  the  plots  are:  the  calculated  temperatures  for  the  furnace-exposed 
surface  of  panel  1,  designated  as  Exp.;  the  temperature  of  the  furnace,  Tstd’  as  determined  from  Eq. 
(29);  and  the  temperatures  calculated  at  locations  corresponding  to  the  locations  BL/Cav.[exp.], 
BL/Cav.[unexp.],  BL/FL[unexp.],  and  UnExp.[under  pads]  ( or  UnExp.[bare]). 

As  can  be  seen  in  Figures  9a-c,  and  as  was  the  case  with  the  1x1  asembly/test  simulations  of  Figures 
6a-c,  when  the  CpuR  and  the  Gj’s  are  0.9,  0.8,  or  0.7,  the  comparisons  between  simulated  and 
measured  temperature  are  very  favorable  (up  to  the  time  when  the  calculation  is  invalidated  because 
of  panel- 1 failure,  when  the  measured  temperature  at  BL/Cav.[exp.]  reaches  approximately  600  °C), 
with  comparisons  improving  with  increasing  emissivity.  Of  the  three  simulations,  the  two 
corresponding  to  6puR  = €;  = 0.9  and  0.8  (Figures  9a  and  9b)  led  to  the  best  comparisons  between 
measured  and  calculated  values,  again  with  little  variation  between  the  results  in  these  two  cases. 
Also,  in  the  0.9  - 0.7  range  of  emissivity,  the  dependence  of  simulated  temperatures  on  emissivity 
is  again  seen  to  be  relatively  weak,  where,  in  contrast  to  this,  a clear  degradation  in  comparison 
between  measurements  and  calculations  occurs  when  the  Cj’s  are  reduced  to  0.5  (see  Figure  9d). 

In  the  1x2  assembly/test,  where  data  from  individual  thermocouples  are  available  from  [1],  it  is 
interesting  to  get  a sense  of  the  relative  variation  of  the  different  thermocouple  measurements  that 
make  up  the  previously-plotted  averages  of  the  measured  temperatures.  Toward  this  end,  all 
measured  temperatures  at  the  generic  locations  BL/Cav.[exp.]  and  BL/Cav.[unexp.]  are  plotted  in 
Figure  10  together  with  the  Figure  9a  calculated  values.  It  is  seen  from  this  plot  that,  for  the  1x2 
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assembly/test  and  throughout  the  test  times  of  interest,  the  characteristic  variations  of  measured 
temperatures  at  ca\'it>’  surface  2 and  at  cavity  surface  3 are  the  same  as  the  characteristic  difference 
between  the  spacially-averaged  temperatures  of  these  two  surfaces.  The  plot  also  reveals  clearly  the 
development  of  the  failure  of  panel  1 as  the  temperature  throughout  its  depth  uniformly  exceeds 
approximately  600  °C. 

Calculation  times  for  all  above  1x2  assembly/test  simulations  were  approximately  270  s. 

Simulations  with  NPTS  = 11:  Figure  6-type  simulations  were  carried  out  with  all  parameters, 
except  for  NPTS,  identical  to  those  used  previously.  Here,  NPTS  was  reduced  to  1 1. 

Calculation  times  for  all  the  new  simulations  were  reduced  from  approximately  270  s to  45  s.  The 
results  of  the  particular  simulation  that  corresponds  to  that  of  Figure  9a  are  presented  in  Figure  11. 
As  can  be  seen  from  Figures  9a  and  1 1,  and  as  was  the  case  in  the  1x1  assembly/test  simulations  the 
calculated  temperatures  for  this  new  reduced-NPTS  calculation  do  not  compare  as  favorably  to 
measured  values  as  do  the  NPTS  = 21  calculation  results.  Again,  use  of  the  relatively-coarse  NPTS 
grid  has  introduced  “wiggles,”  here  even  more  pronounced  then  for  the  1x1  test,  into  the  calculated 
temperamre  histories.  And  again,  at  significant  savings  in  computational  time,  the  new  calculation 
provides  results  which  still  compare  relatively  well  with  the  measured  temperatures. 

In  the  case  of  the  Figure  9b-d-t>pe  simulations,  the  comparison  of  the  NPTS  = 1 1 to  the  NPTS  = 21 
results  were  similar  to  the  Figure  9a/Figure  1 1 result. 

GYPST  Simulations  With  “Nearly”-Discontinuous  Density  and  Thermal  Conductivity 
Representations 

As  discussed  and  as  presented  in  APPENDIX  A and  in  Figures  3-5,  all  the  above  GYPST 
simulations  use  gypsum  properties  where  Reference- [2] -recommended,  discontinuous,  step  changes 
of  p(T)  and  k(T),  at  T = 80  °C  and  100  °C,  respectively,  are  represented  by  corresponding  changes 
over  a 20  °C  interval. 

It  is  of  interest  to  study  the  effect  on  a GYPST  simulation  of  using  property  representations  that  are 
“more-nearly”  discontinuous,  e.g.,  more  in  keeping  with  the  Reference  [2]  recommendations,  with 
significant  changes  in  properties  being  represented  as  occuring  over  intervals  of  only  a few  °C. 
Toward  this  end.  one  additional  calculation  was  carried  out,  identical  to  that  of  Figure  9a,  except  that 
temperature-dependent  changes  of  p(T)/p(T  = 20  °C)  and  k(T),  centered  at  T = 80  °C  and  100 °C, 
respectively,  and  taking  place  over  a 20  °C  interval,  are  now  represented  as  taking  place  over  a 2 °C 
intern  al  (i.e.,  70  and  90  of  Table  A-2  are  replaced  by  79  and  81,  respectively,  and  the  90  and  1 10  of 
Table  A-4  are  replaced  by  99  and  101,  respectively). 

The  result  of  the  new  simulation  is  presented  in  Figure  12,  which  is  to  be  compared  to  Figure  9a. 
As  can  be  seen,  except  for  the  calculated  temperature  at  the  thermocouple  location  BL/Cav.[unexp.] 
(i.e.,  at  the  cavity-exposed  surface  of  panel  2),  and  then  only  when  temperatures  there  are  less  than 
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400  °C,  differences  in  the  two  results  are  hardly  noticeable. 

Compared  to  the  calculation  time  of  the  Figure  9a  simulation,  calculation  time  for  that  of  Figure  12 
was  unchanged  at  approximately  270  s. 

The  above  comparisons  indicate  that,  computationally,  GYPST  simulations  are  robust  when  dealing 
with  “nearly”-discontinuous  p(T),  Cp(T),  and  k(T)  representations.  However,  they  also  indicate  that 
some  care  in  representing  the  real  “smoothed”  nature  of  the  temperature  dependence  of  these 
properties  (i.e.,  significant  property  changes  over  several  °C  V5  possibly  unrealistic  “nearly”- 
discontinuous  step  changes)  leads  to  greater  confidence  in  the  uniform  accuracy  of  a simulation. 


SUMMARY  AND  CONCLUSIONS 

This  work  developed  a methodology  and  an  associated  FORTRAN  subroutine,  GYPST,  for  solving 
an  initial-value/boundary-value  problem  to  simulate  the  thermal  response  of  fire-environment- 
exposed  wall  systems  constructed  of  two  gypsum  panels  of  arbitrary  thickness  mounted  one  on  each 
side  of  an  array  of  steel  studs.  GYPST  calculates  time-dependent  temperatures  through  the  thickness 
of  the  panels  at  equidistant  sets  of  node  points. 

The  wall  model  was  verified  by  using  GYPST,  in  its  stand-alone  mode,  to  simulate  the  ASTM  El  19 
fumace-fire  test  performance  of  two  different  wall  system  designs,  where  results  of  furnace  tests  on 
the  two  designs  are  reported  in  [1]  and  [2].  Calculated  temperature  response  of  the  gypsum  panel 
surfaces  for  both  of  the  wall  designs  compared  favorable  with  measured  temperatures  when  GYPST - 
user-specified  panel-surface  and  furnace  emissivities  were  specified  in  the  range  0.7  - 0.9.  Excellent 
results  were  achieved  when  all  emissivities  were  chosen  to  be  0.8  or  0.9,  with  little  variation  in 
results  between  these  two  cases. 

The  ASTM  El  19  simulations  showed  that  very  favorable  comparisons  between  calculation  and 
experiment,  obtained  with  20  or  21  node  points  per  panel  (for  the  two  wall  design  systems),  were 
degraded,  but  were  still  useful,  when  the  number  of  node  points  were  reduced  to  10  or  11.  This 
halving  of  the  number  of  node  points  reduced  computation  times  to  a sixth  of  their  original  values, 
from  270  s to  45  s for  one  of  the  wall  design  systems,  and  from  160  s to  25  s in  the  case  of  the  other. 

In  GYPST,  temperature-dependent  properties  for  the  gypsum  panels  are  determined  by  linear 
interpolation  between  user-supplied  (property,  temperature)  data  pairs.  Data  pairs  used  for  all 
GYPST  simulations  discussed  in  this  work  were  based  on  measured  properties  of  the  gypsum 
material  used  in  the  tested  wall  systems  of  [1]  and  [2]. 

The  product  of  the  above  is  an  experimentally-verified,  modular,  thermal-wall-model  algorithm  and 
associated  computer  subroutine,  GYPST,  that  can  be  integrated  into  zone-type  compartment  fire 
models  and  that  can  also  be  used  for  a variety  of  “stand-alone”  simulations,  including  simulations 
to  determine  the  ASTM  El  19  fumace-fire  test  performance  of  some  wall  systems. 
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NOMENCLATURE 


F,  for  furnace-exposed  surface  1 

F. 

function  for  determining  q',',  Eq.  (2) 

f(GrPr) 

Eq.(12) 

Gr 

Grashoff  number,  Eq.  (12) 

a 

acceleration  of  gravity 

H 

height  of  the  wall  system 

^AMB 

heat  transfer  coefficient,  surface-to-ambient,  Eq.  (5) 

FpL-R 

heat  transfer  coefficient,  surface-to-fumace,  Eq.  (4) 

^EFF 

effective  k of  the  air  gap,  Eq.  (12) 

Lj,  L2, 

Eq.(l) 

Pr 

Prandtl  number  of  air 

q” 

net  rate  of  heat  transfer  to  surface  I 

qi.CON’  qS.CON 

concective  components  of  q^,  q's 

Q-Z’Rad’  Q3.RAD  radiative  components  of  q'l 


T 

temperature 

Fa.MB 

ambient  temperature 

yd)  y'2) 

T in  the  right,  left  panel,  see  Eigure  (2). 

E.ajr 

characteristic  T of  the  air  in  the  gap,  Eq.  (10) 

Eimt 

T(x,  t = t,) 

T, 

T of  surface  i 

TstD 

T(t)  for  a standard  fire 

20 


t 


time 


X transfromed  x,  Eq.  (21) 

X distance,  see  Figure.  (2) 

Xj  X of  surface  i 

£puj^  effective  emissivity  of  furnace  walls  and  gases 

£;  emissivity  of  surface  I 

Ki,  K2  Eqs.  (21),  (22) 

Eqs.  (21),  (22) 

V kinematic  viscosity  of  the  air 

o Stephan-Bolzmann  constant,  Eq.  (4). 


21 


REFERENCES 


[ 1 ] Sultan,  M.  A.,  et  al.  Temperature  Measurements  in  Full-Scale  Insulated  and  Non-Insulated 
(1x2)  Gypsum  Board  Protected  Wall  Assemblies  with  Steel  Studs,  Internal  Rpt.  675, 
National  Research  Council  Canada  (NR)  Institute  for  Research  in  Construction,  1994. 

[2]  Sultan,  M.A.,  A Model  for  Predicting  Heat  Transfer  Through  Non-Insulated  Unloaded  Steel 
Stud  Gypsum  Board  Wall  Assemblies  Exposed  to  Fire,  Fire  Technology,  Vol.  32,  No.  3, 
pp.239-259,  1996. 

[3]  Standard  Methods  of  Fire  Tests  of  Building  and  Construction  Materials,  ASTM  El  19, 
ASTM  Committee  E-5,  Subcommittee  E05.1 1,  American  Society  for  Testing  and  Materials 
(ASTM),  Philadelphia,  August  1988. 

[4]  Fire  Resistance  Tests  - Elements  of  Building  Construction,  ISO-834,  International 
Organization  for  Standardization,  Switzerland,  1992. 

[5]  Eurocode  on  Actions  on  Structures,  Commission  of  the  European  Communities,  Draft,  p.  9, 
April  1990. 

[6]  Milke,  J.,  Analytic  Methods  for  Determining  Fire  Resistance  of  Steel  Members,  Section 
3/Chapter  6,  SFPE  Handbook  of  Fire  Protection  Engineering,  DiNenno,  P.J.  et  al,  Eds.,  1st 
Edition,  Society  of  Fire  Protection  Engineers,  Boston,  p.  3-100,  1988. 

[7]  Baumeister,  T.  And  Marks,  L.S.,  Eds.,  Standard  Handbook  for  Mechanical  Engineers,  7th 
Edition,  McGraw  Hill,  p.  4-103,  1958. 

[8]  Eckert,  E.R.G.,  and  Drake,  R.M.,  Heat  and  Mass  Transfer,  2nd  Edition,  McGraw  Hill,  pp. 
331-332,  1959. 

[9]  Kraussold,  H.,  et  al,  Forschung  a.d.  Gebiete  d.  Ingenieurw,  5,  pp.  186-191,  1934. 

[10]  Jacob,  M.,  Heat  Transfer.  Vol.  1,  John  Wiley,  p.  542,  1949. 

[11]  Mehaffey,  J.R.,  Cuerrier,  P.,  and  Carisse,  G.,  A Model  for  Predicting  Heat  Transfer  Through 
Gypsum-BoardAVood-Stud  Walls  Exposed  to  Fire,  Fire  and  Materials,  18,  pp.  297-305, 
1994. 

[12]  Hurst,  J.P.,  and  Ahmed,  G.N.,  Modeling  the  Thermal  Response  of  Gypsum  Wallboard  and 
Stud  Assemblies  Subject  to  Standard  Fire  Testing,  Proceedings  of  the  1995  International 
Conference  on  Fire  Research  and  Engineering,  Orlando,  D.P.  Lund,  Ed.,  Society  for  Fire 
Protection  Engineers,  Boston,  pp.  568-573,  1995. 


22 


[ 1 3]  Mehaffey,  J.R.  and  Tekeda,  H,  Predicting  Heat  Transfer  Through  Wood-Stud  Walls  Exposed 
to  Fire,  Proceedings  of  the  1995  International  Conference  on  Fire  Research  and  Engineering. 

Orlando,  D.P.  Lund,  Ed.,  Society  for  Fire  Protection  Engineers,  Boston,  pp.  557-562,  1995. 

[14]  Gerlich,  J.T.,  Collier,  P.C.R.,  and  Buchanan,  A.H.,  Design  of  Light  Steel-framed  Walls  for 
Fire  Resistance,  Fire  and  Materials,  20,  pp. 79-96,  1996. 

[15]  Collier,  P.C.R.,  A Model  for  Predicting  the  Fire-Resisting  Performance  of  Small-Scale 
Cavity  Walls  in  Realistic  Fires,  Fire  Technology,  pp.  120-136,  2nd  Quarter,  1996. 

[16]  Thomas,  G.C.,  Buchanan,  A.H.,  Carr,  A.J.,  Fleishmann,  C.M.,  and  Moss,  P.J.,  Light  Timber- 
Framed  Walls  Exposed  to  Compartment  Fires,  Jl.  Fire  Protection  Engineering,  7(1),  pp.25- 
35,  1995. 

[17]  Hyman,  J.M.,  MOL  ID,  Rev.  1,  Los  Alamos  Scientific  Laboratory,  Los  Alamos,  Oct.  1978. 

[18]  Hindmarsh,  A.C.,  GEARB:  Solution  of  Ordinary  Differential  Equations  Having  Banded 
Jacobian,  University  of  California  Report  UCID-30059,  March,  1975. 


23 


APPENDIX  A: 


THE  THERMAL  PROPERTIES  REPRESENTATIONS  USED  IN  THE 
WALL  MODEL 


Linear  Interpolation  Between  Specified  Pairs  of  Property,  Temperature  Data  Points 

This  Appendix  presents  the  analytic  representations  used  in  the  wall  model  for  the  properties 
P(T)/p(Tr£p),  Cp(T),  and  k(T)  of  the  gypsum  panels.  All  these  are  all  defined  by  an  algorithm, 
incorporated  in  the  GYPST  software,  that  involves  linear  interpolation  in  temperature  between 
specified  [property,  temperature]  data  pairs  and  extrapolation  at  constant  property  beyond  the 
extreme  temperatures  of  these  data  pairs.  For  any  particular  property,  the  only  requirement  for  a 
specified  set  of  data  pairs  is  that  it  be  consistant  with  property  values  that  are  continous  in 
temperature  (i.e.,  step  discontinuities  are  not  allowed). 

While  particular  gypsum  board  products  can  have  a significant  differences  in  ambient-  temperature 
density.  Reference  [11]  suggests  that  for  any  reference  temperature,  Tj^p,  the  dimensionless  density 
function  p(T)/p(TR£p)  is  approximately  invarient  from  one  product  to  another.  Based  on  this  idea, 
the  representation  in  GYPST  for  p(T)  of  a particular  product  is  obtained  by  the  specified  value 
p(T  = 20  °C)  of  the  product  and  by  knowledge  of  the  assumed-universal  density  function 
p(T)/p(T  =20°C).  Thus,  if  a representation  of  p(T)/p(T  = 20  °C)  is  obtained  from  one  set  of 
measurements  on  one  particular  gypsum  product,  there  is  a basis  for  using  it  to  estimate 
p(T)/p(T  = 20  °C)  for  a wide  range  of  other  products. 

Properties  Representations  Used  in  GYPST  Simulations  of  the  Tests  of  References  [1]  and  [2] 

Along  with  the  measured  value 


p(20  °C)  = 698  kg/m^ 


(Al) 


the  following  data  pairs  are  provided  by  Reference  [2],  to  represent,  via  linear  interpolation, 
measured  values  of  p(T)/p(TREp),  Cp(T),  and  k(T)  for  the  gypsum  panels  used  in  the  tested  wall 
systems  of  [1]  and  [2]: 


Density: 


p(T)/p(T  = 20  °C) 


T [°C] 


(A-2) 


1. 

0.825 


70 

90 
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Specific  Heat: 


Cp  [J/(kgK)]  T [°C] 


1500 

20 

1842 

78 

2769 

85 

5861 

97 

18479 

124 

2006 

139 

1001 

148 

714 

373 

715 

430 

571 

571 

618 

609 

3000 

662 

3070 

670 

571 

685 

(A-3) 


Thermal  Conductivity: 


k [W/(m  K)]  T [°C] 


0.25 

90 

0.12 

no 

0.12 

370 

0.27 

800 

1.83 

2000 

(A-4) 


Comment:  Note  that  in  Reference  [1],  p(T)/p(T  = 20  °C)  and  k(T)  are  represented  as  being 
discontinuous  (i.e.,  step  changes)  at  T = 80  °C  and  100  °C,  respectively.  These  approximations  are 
represented  here  as  changes  that  takes  place  continuously  over  20  °C  temperature  interv'als  centered 
at  80°C  [see  (A-2)]  and  100  °C  [see  (A-4)],  respectively.  Other  similar,  but  less  obvious 
discontinuities  also  exist  in  the  reference-[2]  representation  of  Cp(T)  and  k(T).  These  discontinuities 
are  removed  in  the  present  representations  of  these  functions,  i.e.,  linear  interpolation  between  the 
data  pairs  of  Eqs.  (A-3)  and  (A-4). 

Plots  of  p(T)/p(T  = 20  °C),  Cp(T),  and  k(T)  from  Eqs.  (A-l)-(A-4)  and  the  interpolation/extrapolation 
algorithm  of  GYPST  are  presented  in  Figures  3-5,  respectively. 
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APPENDIX  B:  A DESCRIPTION  OF  THE  GENERAL  USAGE  OF  GYPST  (FROM 

THE  GTPSr  SOURCE  CODE) 


SUBROUTINE  GYPSTfTEMPI ,TEMP2,IXM,XM1  ,XM2,NPTS,IBCTYP,EPS1  ,EPS2, 

1 EPS3,EPS4,EPSFUR,EPSCAL,TINT,TLAST,TMP1  IN,TMP4IN,TMP1  FN 

2 ,TMP4FN,TAMB,Q1  IN,Q1  FN,Q4IN,Q4FN,AL1  ,AL2,ALGAP,DENAMB 

3 ,MRHOD,RHODDAT,MCP,CPDAT,MCOND,CONDDAT) 

0 «*«««*««*«««**•************«*««****«■**«  ***««***«***«-*«*«*******«**•***** 

C THIS  PROGRAM  INTEGRATES  THE  INITIA17BOUNDARY  VALUE  PROBLEM,  BETWEEN  TIME=TINT  AND 

C TIME=TLAST,  FOR  THE  TEMPERATURE  OF  A PAIR  OF  GYPSUM  PANELS  SEPARATED  BY  AN  AIR  GAP 

C AND  HEATED  ON  EITHER  EXPOSED  SURFACE.  INPUT  INCLUDES  THE  TEMPERATURE  DISTRIBUTION 

C AT  TIME=TINT  AND  AT  NPTS  EQUIDISTANT  POINTS  IN  EACH  OF  THE  TWO  PANEL,  INCLUDING  A POINT 

C AT  EACH  PANEL  SURFACE.  OUTPUT  INCLUDES  THE  NEW  TEMPERATURE  DISTRIBUTION  AT  THE 

C SAME  POINTS  AND,  IN  THE  CASE  OF  BOUNDARY  CONDITIONS  WHERE  OUTER  SURFACE  HEAT 

C TRANSFER  RATES  ARE  NOT  EXPLICITLY  SPECIFIED,  THE  NEW  VALUES  FOR  THESE  HEAT  TRANSFER 

C RATES. 

C 

C DEFINITIONS: 

C 

C AL1 

C THICKNESS  OF  LEFT  PANEL  [M] 

C AL2 

C THICKNESS  OF  RIGHT  PANEL  [M] 

C ALGAP 

C THICKNESS  OF  AIR  GAP  [M] 

C CONDDATA(1  ,N)=CONDDAT(1  ,N),  N=1  ,NCOND=MCOND 

C COND  [W/(M*K)]  OF  NTH  (COND,TEMP)  DATA  PAIR  USED  TO  CALCULATE  COND(TEMP)  IN 

C SUBROUTINE  CONDGYP 

C CONDDATA(2,N)=CONDDAT(2,N),  N=1,NCOND=MCOND 

C TEMP  [C]  OF  NTH  (COND, TEMP  DATA  PAIR  USED  TO  CALCULATE  COND(TEMP)  IN 

C SUBROUTINE  CONDGYP 

C CPDATA(1  ,N)=CPDAT(1  ,N),  N=1 , NCP=MCP 

C CP  [J/{KG*K)]  OF  NTH  (CP,TEMP)  DATA  PAIR  USED  TO  CALCULATE  CP(TEMP)  IN  SUBROUTINE 

C DSPGYP 

C CPDATA(2,N)=CPDAT(2,N),  N=1,  NCP=MCP 

C TEMP  [C]  OF  NTH  (CP, TEMP)  DATA  PAIR  USED  TO  CALCULATE  CP(TEMP)  IN  SUBROUTINE 

C DSPGYP 

C DENAM 

C RHO  OR  DENSITY  OF  GYPSUM  AT  20  C [KG/M**3] 

C EPS1,EPS2,EPS3,AND  EPS4 

C THE  EMISSIVITIES  OF  GYPSUM  SURFACES  1 , 2,  3,  AND  4,  RESPECTIVELY. 

C EPSFUR 

C EFFECTIVE  EMISSIVITY  OF  THE  FURNACE,  FOR  A STANDARD  FIRE  EXPOSURE  OF  SURFACE 

C 1 OR  OF  SURFACES  1 AND  4,  OR  THE  EFFECTIVE  EMISSIVITY  OF  THE  FIRE/FACING- 

C ROOM-SURFACES  FOR  ROOM  FIRE  EXPOSURES  OF  SURFACE  1 OR  SURFACES  1 AND  4. 

C EPSCAL 

C THE  SUBROUTINE  USED  TO  CARRY  OUT  THE  INTEGRATION  IS  MOL1 D AND  ITS  ASSOCIATED 

C SUBROUTINES.  IN  THIS,  EPSCAL  IS  THE  TIME  INTEGRATOR  LOCAL  ERROR  TOLERANCE. 

C PARAMETER  ESTIMATES  OF  THE  TIME  STEP  ERROR  AT  THE  SPATIAL  GRID  POINTS  FOR  ALL 

C GRID  POINTS  FOR  ALL  THE  PDES  IS  KEPT  LESS  THAN  EPSCAL  IN  THE  ROOT-MEAN-SQUARE 

C (RMS)  NORM.  EPSCALIS  ALSO  USED  TO  CHOOSE  THE  INITIAL  STEP. 

C FIRE1DAT(1,N),  N=1,  NFIRE1 

C TEMPERATURE  [K]  OF  NTH  DATA  PAIR  (TEMPERATURE,TIME)  USED  IN  SUBROUTINE  FRTMP1 

C TO  CALCULATE  TMPFR1(TIME),  THE  TEMPERATURE  VS  TIME  OF  A FIRE  ENVIRONMENT  TO 

C WHICH  SURFACE  1 IS  EXPOSED  WHEN  IBCTYP  = 7 OR  8;  NFIRE1  =<  1000. 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


FIRE1DAT(2,N),  N=1,  NFIRE1 

TIME  [S]  OF  NTH  DATA  PAIR  (TEMPERATURE, TIME)  USED  IN  SUBROUTINE  FRTMP1  TO 
CALCULATE  TMPFR1  (TIME),  THE  TEMPERATURE  VS  TIME  OF  A FIRE  ENVIRONMENT  TO 
WHICH  SURFACE  1 IS  EXPOSED  WHEN  IBCTYP  = 7 OR  8;  NFIRE1  =<  1000. 

FIRE4DAT(1  ,N),  N=1 , NFIRE4 

TEMPERATURE  [K]  OF  NTH  DATA  PAIR  (TEMPERATURE, TIME)  USED  IN  SUBROUTINE  FRTMP4 
TO  CALCULATE  TMPFR4(TIME),  THE  TEMPERATURE  VS  TIME  OF  A FIRE  ENVIRONMENT  TO 
WHICH  SURFACE  4 IS  EXPOSED  WHEN  IBCTYP  = 8;  NFIRE4  =<  1000. 

FIRE4DAT(2,N),  N=1,  NFIRE4 

TIME  [S]  OF  NTH  DATA  PAIR  (TEMPERATURE,TIME)  USED  IN  SUBROUTINE  FRTMP4  TO 
CALCULATE  TMPFR4(TIME),  THE  TEMPERATURE  VS  TIME  OF  A FIRE  ENVIRONMENT  TO 
WHICH  SURFACE  4 IS  EXPOSED  WHEN  IBCTYP  = 8;  NFIRE4  =<  1000. 

IBCTYP 


IXM 


TYPE  OF  BOUNDARY  CONDITIONS: 

IBCTYP  = 1:  [01  IN,  01 FN]  AND  [04IN,  04FN]  SPECIFIED;  01  SPECIFIED  AS  A LINEAR 

INTERPOLATION  IN  TIME  BETWEEN  01  IN,  AT  TIME  TINT,  AND  01  FN,  AT  TIME 
TLAST;  04  SPECIFIED  SIMILARLY  BY  04IN  AND  04FN. 

IBCTYP  = 2:  [01  IN,  01  FN]  AND  [TMP4IN,  TMP4FN]  SPECIFIED;  01  SPECIFIED  AS  LINEAR 

INTERPOLATION  IN  TIME  BETWEEN  01  IN,  AT  TIME  TINT,  AND  01  FN,  AT  TIME 
TLAST;  TEMPERATURE  OF  SURFACE  4 SPECIFIED  SIMILARLY  BETWEEN 
TMP4IN  AND  TMP4FN. 


IBCTYP  = 3: 


IBCTYP  = 4: 


IBCTYP  = 5: 


IBCTYP  = 6: 


IBCTYP  = 7: 


IBCTYP  = 8: 


IBCTYP  = 9: 


[TMP1IN,TMP1FN]  AND  [04IN,  04FN]  SPECIFIED;  TEMPERATURE  OF 
SURFACE  4 SPECIFIED  AS  A LINEAR  INTERPOLATION  IN  TIME  BETWEEN 
TMP1IN,  ATTIME  TINT,  AND  TMP1FN,  AT  TIME  TLAST;  04  SPECIFIED 
SIMILARLY  BY  04IN  AND  04FN. 

[TMP1IN,  TMP1FN]  AND  [TMP4IN,  TMP4FN]  SPECIFIED;  TEMPERATURE  OF 
SURFACE  1 SPECIFIED  AS  A LINEAR  INTERPOLATION  IN  TIME  BETWEEN 
TMP1  IN,  AT  TIME  TINT,  AND  TMP1  FN,  AT  TIME  TLAST;  TEMPERATURE  OF 
SURFACE  4 SPECIFIED  SIMILARLY  BETWEEN  TMP4IN  AND  TMP4FN. 

01  SPECIFIED  AS  HEAT  TRANSFER  FROM  AN  ASTM  El  19  STANDARD-FIRE 
FURNACE  ENVIRONMENT  (TEMPERATURE-TIME  ENVIRONMENT  IS 
TMPSTD(TIME))  AND  04  SPECIFIED  AS  HEAT  TRANSFER  FROM  AN  AMBIENT 
TEMPERATURE  (TAMB)  LABORATORY  ENVIRONMENT. 

01  SPECIFIED  AS  HEAT  TRANSFER  FROM  AN  ISO  834  STANDARD-FIRE 
FURNACE  ENVIRONMENT  (TEMPERATURE-TIME  ENVIRONMENT  IS 
TMPSTI(TIME))  AND  04  SPECIFIED  AS  HEAT  TRANSFER  FROM  AN  AMBIENT 
TEMPERATURE  (TAMB)  LABORATORY  ENVIRONMENT. 

01  IS  THE  HEAT  TRANSFER  FROM  A SPECIFIED  ROOM-FIRE  ENVIRONMENT 
(TEMPERATURE-TIME  OF  THE  FIRE  ENVIRONMENT  IS  TMPFRI(TIME),  CAL- 
CULATED IN  THE  SUBROUTINE  FRTMP1  BY  INTERPOLATION/EXTRAPO- 
LATION OF  DATA  IN  INPUT  COMMON  BLOCK  FIRE1BLK)  AND  04  IS 
SPECIFIED  AS  HEAT  TRANSFER  FROM  AN  AMBIENT-TEMPERATURE  (TAMB) 
LABORATORY  ENVIRONMENT. 

01  IS  THE  HEAT  TRANSFER  FROM  ONE  SPECIFIED  ROOM-FIRE 
ENVIRONMENT  (TEMPERATURE-TIME  OF  THE  FIRE  ENVIRONMENT  IS 
TMPFR1(TIME),  CALCULATED  IN  THE  SUBROUTINE  FRTMP1  BY  INTERPO- 
LATION/EXTRAPOLATION OF  DATA  IN  INPUT  COMMON  BLOCK  FIRE1BLK) 
AND  04  THE  HEAT  TRANSFER  FROM  ANOTHER  SPECIFIED  ROOM-FIRE 
ENVIRONMENT  (TEMPERATURE-TIME  OF  THE  FIRE  ENVIRONMENT  IS 
TMPFR4(T),  CALCULATED  IN  THE  SUBROUTINE  TRTMP4  BY  INTERPOLATION/ 
EXTRAPOLATION  OF  DATA  IN  INPUT  COMMON  BLOCK  FIRE4BLK). 
[TMP1IN,TMP1FN]  SPECIFIED;  THE  TEMPERATURE  OF  SURFACE  1 IS 
SPECIFIED  AS  AS  A LINEAR  INTERPOLATION  IN  TIME  BETWEEN  TMP1IN,  AT 
TIME  TINT,  AND  TMP1  FN,  AT  TIME  TLAST;  04  SIMULATED  TO  BE  THE  HEAT 
TRANSFER  TO  SURFACE  4 FROM  AN  AMBIENT  TEMPERATURE  (TAMB) 
ENVIRONMENT. 
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DETERMINES  WHETHER  XMI(NPTS),  XM2(NPTS)  IS  INPUT  OR  OUTPUT: 

IXM  = 1 IF  XM1,  XM2  ARE  TO  BE  CALCULATED  AND  THEN  PROVIDED  AS  OUTPUT; 
IXM  = 2 IF  XM1  AND  XM2  WERE  CALCULATED  PREVIOUSLY  AND  ARE  BEING 
PROVIDED  AS  INPUT. 

NPTS 

NUMBER  OF  POINTS  IN  PANELS  WHERE  TEMPERATURE  IS  COMPUTED. 

NCOND  = MCOND 

NUMBER  OF  (CONDJEMP)  PAIRS  USED  TO  CALCULATE  COND  VS 
TEMPERATURE  FOR  THE  GYPSUM  PANELS,  =<  100. 

NCP=MCP 

NUMBER  OF  (CP, TEMP)  PAIRS  USED  TO  CALCULATE  CP  VS  TEMPERATURE  FOR  THE 
GYPSUM  PANELS,  =<  100. 

NFIRE1 

NUMBER  OF  (TEMP,  TIME)  PAIRS  USED  TO  CALCULATE  THE  TEMPERATURE-TIME  CURVE 
OF  THE  FIRE  TMPFR1(T)  CALCULATED  IN  THE  SUBROUTINE  FRTMP1,  =<  1000. 

NFIRE4 

NUMBER  OF  (TEMP,  TIME)  PAIRS  USED  TO  CALCULATE  THE  TEMPERATURE-TIME  CURVE 
OF  THE  FIRE  TMPFR4(T)  CALCULATED  IN  THE  SUBROUTINE  FRTMP4,  =<  1000. 

NRHOD  = MRHOD 

NUMBER  OF  (RHO,TEMP)  PAIRS  USED  TO  CALCULATE  RHO  VS  TEMPERATURE  FOR  THE 
GYPSUM  PANELS,  =<100. 

Q1 

RATE  OF  HEAT  TRANSFER  TO  SURFACE  1 [W/M**2]. 

Q4 

RATE  OF  HEAT  TRANSFER  TO  SURFACE  4 [W/M*'‘2]. 

Q1IN,  01  FN 

INITIAL,  FINAL  SPECIFIED  Q1  [W/M**2]. 

Q4IN,  Q4FN 

INITIAL  FINAL  SPECIFIED  Q4  [W/M**2]. 

RHODDATA(1,N)=RHODDAT(1,N),  N=1,  NRHO  = MRHO 

DENSITY  [KG/M**3]  OF  NTH  (DENSITY,TEMP)  DATA  PAIR  USED  TOCALCULATE  RHO(TEMP) 
IN  SUBROUTINE  DSPGYP. 

RHODDATA(2,N)  = RHODDAT(2,N),  N=1,  NRHO=MRHO 

TEMP  [C]  OF  NTH  (DENSITY,TEMP)  DATA  PAIR  USED  TO  CALCULATE  RHO(TEMP)  IN 
SUBROUTINE  DSPGYP. 

SURFACE  1 

OUTER/LEFT  SURFACE  OF  LEFT  PANEL. 

SURFACE  2 

INNER/RIGHT  SURFACE  OF  LEFT  PANEL. 

SURFACE  3 

INNER/LEFT  SURFACE  OF  RIGHT  PANEL. 

SURFACE  4 

OUTER/RIGHT  SURFACE  OF  RIGHT  PANEL. 

TEMPI 

TEMPERATURE  DISTRIBUTION  IN  THE  LEFT  PANEL  (I.E.,  AT  THE  NPTS  POINTS)  AT  TIME  TINT 
[K]. 

TEMP2 

TEMPERATURE  DISTRIBUTION  IN  THE  RIGHT  PANEL  (I.E.,  AT  THE  NPTS  POINTS)  AT  TIME  TINT 
[K]. 

TINT,  TLAST 

INITAL,  FINAL  TIMES  [S]. 

TMP1IN,  TMP1FN 

INITIAL,  FINAL  TEMPERATURES  OF  SURFACE  1 [K]. 

TMP4IN,  TMP4FN 

INITIAL,  FINAL  TEMPERATURES  OF  SURFACE  4 [K]. 

TAMB 

AMBIENT  TEMPERATURE  OF  LABORATORY/OUTSIDE  ENVIRONMENT  [K];  REQUIRED  IF 
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oooooooooooooooooooooooooooooooooooooooooooooooooo 


IBCTYP  = 5,  6,  7,  OR  9. 

X 

DISTANCE  FROM  LEFT  TO  RIGHT  AS  MEASURED  FROM  SURFACE  1. 

XM1(N) 

VALUE  OF  X IN  LEFT  PANEL,  PANEL  1,  AT  ITS  NTH  TEMPERATURE  CALCULATED  POSITION 
FROM  THE  LEFT  [I.E.,  WHERE  TEMPERATURE  = TEMPI (N)  AT  TIME  TINT],  N = 1 , NPTS,  WHERE 
THE  NPTS  ARE  EQUALLY  SPACED  THROUGH  HE  PANEL  THICKNESS  [XM1(1)  = 0.  IS  AT 
SURFACE  1 AND  XMI(NPTS)  = AL1  IS  AT  SURFACE  2]. 

XM2(N) 

VALUE  OF  X IN  RIGHT  PANEL,  PANEL  2,  AT  ITS  NTH  TEMPERATURE-CALCULATED  POSITION 
FROM  THE  LEFT  [I.E.,  WHERE  TEMPERATURE  = TEMP2(N)  AT  TIME  TINT],  N = 1 , NPTS,  WHERE 
THE  NPTS  ARE  EQUALLY  SPACED  THROUGH  THE  PANEL  THICKNESS  [XM2(1)  = AL1  + ALGAP 
IS  AT  SURFACE  3 AND  XM2(NPTS)  = AL1  + ALGAP  + AL2  IS  AT  SURFACE  4]. 

INPUT  (IN  ORDER  OF  CALL  TO  THIS  SUBROUTINE): 

1 ) AT  TIME=TINT:  TEMPI  (N),  TEMP2(N),  N = 1 , NPTS; 

2)  IXM; 

3)  XM1  (N),  XM2(N),  N = 1 , NPTS;  ONLY  REQUIRED  IF  IXM  > 1 ; 

4)  NPTS; 

5)  IBCTYP; 

6)  EPS1,EPS2,EPS3,AND  EPS4; 

7)  EPSFUR;  ONLY  REQUIRED  IF  IBCTYPE=5,  6,  7,  OR  8; 

8)  EPSCAL; 

9)  TINT,  TLAST; 

1 0)  TMP1  IN;  ONLY  REQUIRED  IF  IBCTYP  = 3,  4,  OR  9; 

1 1 ) TMP4IN;  ONLY  REQUIRED  IF  IBCTYP  = 2 OR  4; 

1 2)  TMP1 FN;  ONLY  REQUIRED  IF  IBCTYP  = 3,  4,  OR  9; 

13)  TMP4FN;  ONLY  REQUIRED  IF  IBCTYP  = 2 OR  4; 

14)  TAMB;  ONLY  REQUIRED  IF  IBCTYP  = 5,  6,  7,  OR  9; 

15)  Q1IN,  Q1FN;  ONLY  REQUIRED  IF  IBCTYP  = 1 OR  2; 

1 6)  Q4IN,  Q4FN;  ONLY  REQUIRED  IF  IBCTYP  = 1 OR  3; 

17)  AL1,AL2,  ALGAP; 

18)  DENAMB; 

19)  MRHOD,RHODDAT(2,MRHOD),MCP,CPDAT(2,MCP),MCOND, 

CONDDAT(2,MCOND);  MRHOD  =<  100,  MCP  =<  100,  MCOND  =<  100; 

20)  COMMON/FIRE1/NFIRE1  ,FIRE1  DAT(2,1 000),  ONLY  REQUIRED  IF 
IBCTYP  = 7 OR  8;  NFIRE1  =<1000; 

21)  COMMON/FIRE4/NFIRE4,FIRE4DAT(2,1000),  ONLY  REQUIRED  IF 
IBCTYP  = 8;  NFIRE4  =<  1000. 

OUTPUT: 

1 ) XM1  (N),  XM2(N),  N = 1 , NPTS; 

2)  AT  TIME=TLAST: 

A)  TEMPI (N),  TEMP2(N),  N = 1 , NPTS; 

B)  Q1  FN  = Q1 , THIS  IS  ONLY  COMPUTED  IF  IBCTYP  = 3, 

4,  5,  6,  7,  OR  8; 

C)  Q4FN  = Q4,  THIS  IS  ONLY  COMPUTED  IF  IBCTYP  = 2, 

4,  5,  6,  7,  8,  OR  9. 
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TABLE  1.  Parameters  of  the  1x1  and  1x2  assemblies/tests  and  their  GYPST  simulations. 


1x1  assemblv/test  1x2  assemblv/test 

Assemblv/test  parameters: 

Tamb 

Properties  of  gypsum: 

23  X 26  X 

p(20  °C) 

698  kg/m^ 

p(T)/p(T  = 20  X),  Cp(T),  and  k(T) 

from  APPENDIX  A 

L, 

0.0159  m 0.0127  m 

L. 

0.0159  m 0.0254  m 

L'GAP 

0.0900  m 

Additional  GYPST  input  parameters: 

IBCTYP 

5 

EPSCAL 

10  ^ 10-",  10  ^ 10^ 

10  ^ 10 ' 

Si,  i = 1 - 4 

0.9,  0.8,  0.7,  or  0.5 

0.9,  0.8,  or  0.7 

NPTS 

10  or  20  11  or  21 
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FIRE-EXPOSED  SIDE  IN  TESTS  OF  REFERENCES  [1,  2] 


Exp, 


BL/Cav.[unexp.] 


1X1 
ASSEMBLY: 


r r p ' r/  r 


\BUCav.[exp.]^:^ 


\ V / 

\ 


1X2 
ASSEMBLY: 


r [ r'/r  -.r 


BL/FL[unexp.]'' 


>,UnExp.[bare] 

UnExp.[under  pads] 


Figure  1.  Sketch  of  example  wall  system  designs  (adopted  from  References  [1]  and  [2]); 
locations  of  calculated/measured  temperatures  in  tests  of  [1]  and  [2]. 
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FIRE  2 
OR 

AMBIENT 


Figure  2.  Geometr>  of  the  arbitrary  two-panel  wall  system  and  a depiction  of  the  x=^X 
transformation  of  Eq.  (22j. 
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p/p(T  = 20  °C) 


Figure 


p(T)/p(T  = 20  °C)  representation  of  APPENDIX  A,  Eq.  (A-2)  used  for  the  wall- 
response  simulations. 
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[ J/(KgK)  ] 


Figure  4.  Cp(T)  representation  of  APPENDIX  A,  Eq.  (A-3)  used  for  the  wall-response 
simulations. 
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k [W/(mK)] 


0.6 


0 200  400  600  800  1000 


T[°C] 


Figure  5.  k(T)  representation  of  APPENDIX  A,  Eq.  (A-4)  used  for  the  wall-response 
simulations. 
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1000 


(aj  Simulation  parameters:  NPTS  = 20;  EPSCAL  = 0.00001;  = all  8;  = 0.9 

Figure  6.  1x1  assembly/test:  simulated  and  averaged  measured  temperatures. 

Tst^ ; Exp.  (calc);  BL/Cav.[exp.]  (calc),  VVV(meas.); 

BL/Cav.[unexp.]  (calc),  TTT(meas);  UnExp.  (calc); 

UnExp.[under  pads]  OOO(meas);  UnExp. [bare]  •••. 
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(Oo)l 


(b)  Simulation  parameters:  NPTS  = 20;  EPSCAL  = 0.00001;  = all  8;  = 0.8 

Figure  6.  1x1  assembly/test:  simulated  and  averaged  measured  temperatures. 

(ConPd)  Tstd ; Exp.  (calc);  BL/Cav.[exp.]  (calc),  VVV(meas.); 

BL/Cav.[unexp.]  ..(calc),  TTT(meas); 

UnExp. (calc);  UnExp.[under  pads]  OOO(meas);  UnExp.[bare]  •••. 
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1000 


(cj  Simulation  parameters:  NPTS  = 20;  EPSCAL  = 0.00001;  e^uR  = all  8j  = 0.7 

Figure  6.  1x1  assembly/test:  simulated  and  averaged  measured  temperatures. 

(Cont’d)  Tstd ; Exp.  (calc);  BL/Cav.[exp.]  (calc),  VVV(meas.); 

BL/Cav.[unexp.]  (calc),  TTT(meas); 

UnExp. (calc);  UnExp.[under  pads]  OOO(meas);  UnExp.[bare]  •••. 
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(d)  Simulation  parameters:  NPTS  = 20;  EPSCAL  = 0.00001;  = 0.8,  all  Sj  = 0.5 

Figure  6.  1x1  assembly/test:  simulated  and  averaged  measured  temperatures. 

(Cont’d)  Tstd ; Exp.  (calc);  BL/Cav.[exp.]  (calc),  VVV(meas.); 

BL/Cav.[unexp.]  (calc),  TTT(meas); 

UnExp. (calc);  UnExp.[under  pads]  OOO(meas);  UnExp.[bare]  ###. 
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(Oo)l 


1000 


Figure  7.  Ixlassembly/test:  simulated  [NPTS  = 10;  EPSCAL  = 0.00001;  = all  e 

0.9]  and  averaged  measured  temperatures.  Tg^D > 

Exp.  (calc);  BL/Cav.[exp.]  .(calc),  VVV(meas.); 

BL/Cav.[unexp.]  (calc),  TTT(meas);  UnExp. (calc); 

UnExp. [under  pads]  OOO(meas);  UnExp. [bare]  ###. 
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Figure  8.  Ixlassembly/test:  simulated  [NPTS  = 20;  EPSCAL  = 0.01;  CpuR  = all  = 0.9] 
and  averaged  measured  temperatures.  Tstd ; Exp. (calc); 

BL/Cav.[exp.]  (calc),  VVV(meas.);  BL/Cav.[unexp.] 

(calc),  TTT(meas);  UnExp. (calc);  UnExp.[under  pads] 

OOO(meas);  UnExp.[bare]  •••. 
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{Oo)l 


Figure  9.  1x2  assembly/test:  simulated  and  averaged  measured  temperatures. 

Tstt) ; Exp.  (calc);  BL/Cav.[exp.]  (calc), 

VVVfmeas.);  BL/Cav.[unexp.]  (calc),  TTT(meas); 

UnExp. (calc);  BL/FL[unexp.] (calc),  □□□  (meas); 

UnExp. [under  pads]  OOO(meas);  UnExp. [bare]  ###. 
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(b)  Simulation  parameters:  NPTS  = 21;  EPS  CAL  = 0.00001;  = all  C;  = 0.8 


Figure  9.  1x2  assembly/test;  simulated  and  averaged  measured  temperatures. 

(Cont’d)  Tstd ; Exp.  (calc);  BL/Cav.[exp.]  (calc), 

VVV(meas.);  BL/Cav.[unexp.]  (calc),  TTT(meas); 

UnExp. (calc);  BL/FL[unexp.] (calc),  □□□  (meas); 

UnExp.[under  pads]  OOO(meas);  UnExp.[bare]  •••. 
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(c)  Simulation  parameters:  NPTS  = 21;  EPSCAL  = 0.00001;  = all  Sj  = 0.7 

Figure  9.  1x2  assembly/test:  simulated  and  averaged  measured  temperatures. 

(Cont’d)  Tstt) ; Exp.  (calc);  BL/Cav.[exp.]  .(calc), 

VVV(meas.);  BL/Cav.[unexp.]  (calc),  TTT(meas); 

UnExp. (calc);  BLyFL[unexp.] (calc),  □□□  (meas); 

UnExp. [under  pads]  OOO(meas);  UnExp. [bare]  ••#. 
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(d)  Simulation  parameters:  NPTS  = 21;  EPSCAL  = 0.00001;  ep^R  = 0.8,  all  8;  = 0.5 

Figure  9.  1x2  assembly/test:  simulated  and  averaged  measured  temperatures. 

(Cont’d)  ; Exp.  (calc);  BLyCav.[exp.]  (calc),  VVV(meas.); 

BL/Cav.[unexp.]  (calc),  TTT(meas); 

UnExp.  (calc);  BL/FL[unexp.]  (calc),  □□□  (meas); 

UnExp.[under  pads]  OOO(meas);  UnExp.[bare]  ••#. 
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Figure  10.  1x2  assembly/test:  simulated  [NPTS  = 21;  EPSCAL  = 0.00001;  = all  8 

0.9]  and  individual  measured  temperatures  at  BL/Cav.[exp.]  and 


BLyCav.[unexp.]. 

Tstt) ; Exp.  (calc);  BL/Cav.[exp.]  (calc), 

VVV(meas.);  BL/Cav.[unexp.]  (calc),  TTT(meas); 

BL/FL[unexp.] (calc);  UnExp. (calc). 
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1000 


Figure  11.  1x2  assembly/test:  simulated  [NPTS  = 11;  EPSCAL  = 0.00001;  e^uR  = all  t, 

0.9]  and  averaged  measured  temperatures.  Ts^d > 

Exp.  (calc);  BL/Cav.[exp.]  VVV(meas.); 


BL/Cav.[unexp.]  (calc),  ▼▼T(meas); 

BL/FL[unexp.] (calc),  □□□  (meas);  UnExp. (calc); 


UnExp.[under  pads]  OOO(meas);  UnExp.[bare]  •••. 
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Figure  12.  1x2  assembly/test:  simulated  [with  “nearly”-discontinuous  p(T)  and  k(T); 

NPTS  = 21;  EPSCAL  = 0.00001;  epuR  = all  8^  = 0.9]  and  averaged  measured 
temperatures.  Tsy^ ; Exp.  (calc); 

BL/Cav.[exp.]  (calc),  VVV(meas.); 

BL/Cav.[unexp.]  (calc),  TTT(meas); 

BL/FL[unexp.] (calc),  □□□  (meas);  UnExp. (calc); 

UnExp. [under  pads]  OOO(meas);  UnExp.[bare]  ###. 
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